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Abstract — Daily  high-resolution  Sea- viewing  Wide  Field-of-view 
Sensor  (SeaWiFS)  images  of  the  central  North  Atlantic  Ocean 
(1998-2003)  show  that  temporal  changes  in  the  absorption  coef¬ 
ficient  of  colored  dissolved  organic  matter  (CDOM)  or  “yellow 
substance”  follow  changes  in  phytoplankton  pigment  absorption 
coefficient  in  time.  CDOM  peaks  (between  January  and  March) 
and  troughs  (late  summer  and  fall)  followed  pigment  peaks  and 
troughs  by  approximately  two  and  four  weeks,  respectively.  This 
phase  shift  is  additional  strong  evidence  that  CDOM  in  the  marine 
environment  is  derived  from  phytoplankton  degradation.  The 
common  assumption  of  linear  covariation  between  chlorophyll 
and  CDOM  is  a  simplification  even  in  this  ocean  gyre.  Due  to  the 
temporal  changes  in  CDOM,  chlorophyll  concentration  estimated 
based  on  traditional  remote  sensing  band-ratio  algorithms  may 
be  overestimated  by  about  10%  during  the  spring  bloom  and  un¬ 
derestimated  by  a  similar  10%  during  the  fall.  These  observations 
are  only  possible  through  use  of  synoptic,  precise,  accurate,  and 
frequent  measurements  afforded  by  space-based  sensors  because 
in  situ  technologies  cannot  provide  the  required  sensitivity  or 
synoptic  coverage  to  observe  these  natural  phenomena. 

Index  Terms — Biogeochemical  cycle,  colored  dissolved  organic 
matter  (CDOM),  ocean  color,  phase  shift,  phytoplankton,  remote 
sensing,  sea  surface  reflectance. 


I.  Introduction 

COLORED  dissolved  organic  matter  (CDOM),  also  com¬ 
monly  called  yellow  substance  or  Gelbstoff,  is  the  col¬ 
ored  fraction  of  the  dissolved  organic  matter  (DOM)  pool  in 
the  sea  [11],  The  importance  of  CDOM  in  carbon  cycling  and 
light  regulation  in  the  marine  environment  is  well  recognized 
(see  review  by  [12]  and  references  therein).  CDOM  strongly  ab¬ 
sorbs  ultraviolet  (UV)  and  blue  light,  and  therefore  affects  phy¬ 
toplankton  photosynthesis  and  the  ecosystem  [22],  CDOM  may 
be  degraded  or  bleached  by  sunlight  [23],  a  transformation  that 
releases  carbon  in  the  form  of  C02  and  CO. 

Despite  recent  progress  in  understanding  the  biogeochem¬ 
istry  of  CDOM  in  marine  waters,  the  origin  of  CDOM  in  the 
open  ocean  has  largely  remained  a  mystery.  Some  argue  that 
CDOM  in  the  sea  is  mainly  of  terrestrial  origin  [7],  while  others 
consider  terrigenous  DOM  to  be  only  a  small  fraction  of  total 
marine  DOM  [19],  [27].  A  consideration  of  loss  processes  leads 
to  the  conclusion  that  <50%  of  the  estuarine  supply  of  CDOM 
remains  on  continental  shelves  [10].  Some  suggest  that  marine 
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Fig.  1.  Western  North  Atlantic  portion  of  the  subtropical  gyre,  with  the  study 
location  marked  by  a  cross  (27.5°N  65°W).  Image  is  a  SeaWiFS  chlorophyll 
concentration  map  overlaid  with  the  2000-  and  4000-m  bathymetry  lines. 
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CDOM  is  a  by-product  of  in  situ  algal  cell  degradation  (e.g., 
[4]),  zooplankton  grazing  [21],  or  microbial  activity  [25],  [26]. 
These  conclusions  are  largely  based  on  laboratory  experiments 
or  limited  field  samples,  from  which  it  is  difficult  to  extrapolate 
to  ocean  basin  scales.  I 

Satellite-borne  sensors  that  measure  the  visible  sea  spectral 
reflectance  (VSSR)  (or  ocean  color)  are  particularly  suited  for 
synoptic  (millions  of  square  kilometers  per  satellite  pass)  and 
frequent  (near  daily)  observations  of  important  bio-optical  prop¬ 
erties  of  the  surface  ocean.  Recent  algorithms  seek  to  partition 
the  VSSR  signal  into  various  components,  for  example  to  accu¬ 
rately  derive  absorption  coefficients  of  phytoplankton  pigments 
(aPh)  and  CDOM/detritus  (ocdm)  (e.g.,  [13],  [16],  and  [29]). 
It  is  difficult  to  partition  ocdm  explicitly  into  absorption  coeffi¬ 
cients  of  CDOM  (ag)  and  detritus  («a)  due  to  the  similarity  of 
their  spectral  shapes.  However,  because  o,c]  usually  contributes 
5%  to  20%  or  less  to  acdm  in  the  deep  open  ocean  [25],  ae 
may  be  represented  by  «cdm  in  studies  of  offshore  processes. 
To  date,  however,  the  relationship  between  satellite-derived  aph 
and  ag  has  not  been  well  documented  even  for  the  open  ocean. 

A  hypothesis  of  general  oceanographic  interest  is  whether 
the  bulk  of  the  CDOM  found  in  the  marine  environment  is  a 
by-product  of  algal  biological  degradation  or  not.  A  close  re¬ 
lationship  between  temporal  changes  in  ag  and  aph  would  help 
understand  the  relationship  between  CDOM  and  phytoplankton. 
However,  in  the  past  it  has  been  impossible  to  explore  this  rela¬ 
tionship  in  the  open  ocean  due  to  limitations  of  in  situ  sampling 
frequency,  the  low  sensitivity  of  ship-based  observations,  and 
the  simple  algorithms  used  to  study  changes  in  ocean  color.  New 
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Fig.  2.  Absorption  coefficients  of  phytoplankton  pigment  (aph443,  m“‘)  and  CDOM  («g443,  m-1)  for  the  oligotrophic  surface  ocean  at  27.5°N  65°W  (Fig.  1), 
derived  from  high-resolution  (1-km)  daily  SeaWiFS  data  between  1998  and  2003  with  a  quasi-analytical  algorithm  [16].  Also  shown  are  the  ag443/api,443  ratio, 
as  well  as  the  ratio  between  chlorophyll  concentration  estimated  from  the  band-ratio  OC4v4  algorithm  [Chl(OC4)]  and  api,443. 


high-quality  satellite  observations  of  the  VSSR  and  increasingly 
sophisticated  algorithms  allow  us  to  revisit  this  basic  hypothesis. 

II.  Data  and  Methods 

To  minimize  the  impact  of  lateral  exchanges  and  terrigenous 
inputs  of  colored  organic  and  other  materials  on  our  conclu¬ 
sions,  our  study  was  conducted  in  the  southern  Sargasso  Sea, 
at  27.5°N  65°W.  This  is  in  the  western  portion  of  the  North 
Atlantic  subtropical  gyre,  where  waters  exhibit  homogeneous 
optical  properties  over  horizontal  distances  exceeding  hundreds 
of  kilometers  (Fig.  1 ).  This  region  features  some  of  the  clearest 
water  in  the  world. 

Data  from  the  Sea-viewing  Wide  Field-of-view  Sensor  (Sea¬ 
WiFS)  collected  during  1998-2003  were  processed  with  the 
SeaDAS4.6  software  released  by  NASA’s  Goddard  Space  Flight 
Center.  The  derived  remote  sensing  spectral  reflectance  data 
(jRrs(A),  sr-1),  after  a  rigorous  quality  control  process  (e.g., 
pixels  with  suspicious  flags  were  removed,  and  a  15  x  15  box 
centered  at  the  location  was  used  to  remove  sensor/algorithm 
noise  [14]),  were  then  used  as  input  to  a  quasi-analytical  al¬ 
gorithm  [16]  to  estimate  aph443  and  ag443  (where  443  is  the 
band  center  wavelength  in  nanometers).  Briefly,  total  absorption 
and  backscattering  coefficients  at  555  nm  were  first  estimated 
empirically  and  analytically.  The  coefficients  at  other  wave¬ 
lengths  were  then  derived  using  an  empirical  estimate  of  particle 
backscattering  spectral  shap'e.  aph443  and  ag443  were  analyt¬ 
ically  derived  from  the  total  absorption  coefficients  at  412  and 
443  nm.  In  this  last  step,  a  spetific  CDOM  absorption  spectral 


slope,  S ,  was  assumed  to  be  constant  and  equal  to  0.02nm_1. 
A  sensitivity  analysis  was  conducted  to  test  this  assumption, 
and  the  results  are  discussed  below.  Chlorophyll-n  concentra¬ 
tion  (Chi)  was  derived  using  the  OC4v4  algorithm  [28]. 

III.  Results  and  Discussion 

Fig.  2  shows  the  SeaWiFS  time  series  of  ag443,  aph443, 
ag443/api,443,  and  the  Chl/aph443  ratio  for  our  Sargasso  Sea 
location.  A  seasonal  cycle  is  evident  in  all  parameters,  with 
winter  maxima  and  late  summer  minima.  There  is  little  interan¬ 
nual  variation,  suggesting  that  the  oligotrophic  ocean  gyre  may 
be  used  to  monitor  long-term  stability  of  future  satellite  sensors. 

The  high  temporal  resolution  and  the  length  of  the  series 
reveal  a  lag  in  ag443  relative  to  aph443.  The  lag  is  more 
manifest  when  all  data  of  1998-2003  are  averaged  to  produce 
a  weekly  climatology  [Fig.  3(a)],  Correlation  analysis  based 
on  the  weekly  climatology  shows  that  ag443  maxima  lagged 
aph443  maxima  by  about  two  weeks  during  winter-spring 
(week  49-19).  The  minimum  ag443  lags  aph443  by  about  4-5 
weeks  during  summer-fall  (week  26-48)  (Fig.  4).  The  lag, 
or  phase  shift,  between  ag443  and  aph443  is  more  manifest 
when  summer-fall  (week  26-48)  lows  and  winter-spring  (week 
49-19)  highs  are  treated  separately,  because  of  the  various 
processes  that  drive  these  temporal  changes  (see  below).  The 
horizontal  lines  in  Fig.  4  correspond  to  the  95%  significance 
leyels  calculated  using  a  “phase  randomization”  method  [8]. 
Correlation  coefficients  above  these  lines  indicate  that  the 
two  time-series  are  significantly  correlated.  The  same  analysis 
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Fig.  3.  (a)  Weekly  climatology  (1998-2003)  of aph443  (per  meter)  and  ng443 
(per  meter)  for  the  oligotrophic  surface  ocean  at  27.5°N  65°W  (Fig.  1).  (b) 
Average  wind  speed  for  the  corresponding  weeks,  obtained  from  the  National 
Center  for  Environmental  Prediction.  Error  bars  indicate  standard  deviation.  In 
order  to  show  the  peaks  and  troughs  more  clearly,  the  week  starts  from  14  to  52, 
and  then  continues  to  week  1. 


Weeks  shifted 

Fig.  4.  Correlation  coefficients  between  aph443  (per  meter)  and  ng443  (per 
meter)  from  the  weekly  climatology  data  [Fig.  3(a)].  Three  analyses  were  per¬ 
formed  for  the  entire  dataset  (week  1-52),  winter-spring  (week  49-19),  and 
summer-fall  (week  26-48),  respectively.  Negative  number  means  that  ag443 
is  behind  api,443.  The  horizontal  lines  are  the  corresponding  95%  significance 
level  for  the  correlation  (see  text  for  more  details). 

could  not  be  applied  to  the  original  daily  data  because  they  were 
unequally  spaced  in  time,  and  interpolation  in  temporal  space 
would  create  additional  uncertainties  due  to  residual  noise  in 
the  daily  data.  Although  the  correlation  coefficient  from  the 
weekly  climatology  would  be  different  than  from  the  daily  data, 
our  focus  is  on  the  relative  changes  versus  time.  Therefore  the 
results  should  be  valid. 

The  phase  shift  suggests  that  the  surface  CDOM  peak  is 
likely  the  result  of  phytoplankton  growth  and  degradation.  This 
may  help  explain  the  findings  of  Nelson  et  al.  [25]  and  Siegel 
and  Michaels  [30]  where  less-frequent  in  situ  data  from  the  U.S. 
Joint  Global  Ocean  Flux  Study’s  Bermuda  Atlantic  Time-Series 
Study  (BATS)  suggested  that  nonliving  materials  (dissolved  and 
particulate)  did  not  appear  to  covary  with  phytoplankton  over 
the  course  of  about  two  years  in  the  upper  60  m.  Clearly,  the 
lag  would  reduce  the  linear  correlation  between  these  variables 
measured  at  the  same  time.  Also,  the  satellite  integrates  the 
signal  over  the  first  optical  depth  in  estimating  the  absorption 


coefficient,  or  about  the  upper  40-50  m  of  the  water  column  at 
this  clear  water  site,  while  a  laboratory  instrument  is  typically 
restricted  to  10-cm  pathlength  sample  cells  to  measure  ag. 
Therefore,  the  satellite  sensor  is  more  sensitive  than  a  laboratory 
instrument,  allowing  for  better  resolution  of  synoptic  changes. 

The  phase  shift  between  ag443  and  aph443  was  manifested 
in  the  seasonality  of  the  ag443/aph443  ratio,  which  suggests 
that  after  a  bloom  CDOM  is  destroyed  at  a  lower  rate  than 
phytoplankton  degradation.  Similarly,  while  the  increase  and 
decrease  in  oph  are  relatively  symmetrical  around  a  bloom, 
ag  buildup  and  decay  are  not.  CDOM  increases  rapidly  and 
peaks  approximately  in  February-March  immediately  after  the 
bloom,  but  then  decays  more  slowly  through  September  and 
October. 

Can  the  similar  aph  and  ag  patterns  be  simply  a  matter  of 
coincidence  with  no  dependence  between  each  other?  Deeper 
wind-driven  and  convective  mixing  [Fig.  3(b)]  in  the  region 
during  winter  brings  deep  nutrients  to  the  surface,  and  the 
availability  of  light  stimulates  higher  phytoplankton  growth 
following  the  classical  Sverdrup  [32]  model.  We  may  speculate 
that  mixing  may,  at  the  same  time,  introduce  deep  unbleached 
CDOM  to  surface  waters.  However,  if  this  is  a  dominant 
process,  we  would  expect  simultaneous  increases  in  both  aph 
and  ag  or  even  a  slight  lead  in  ag  in  the  winter,  since  it  may  take 
time  for  phytoplankton  to  develop  the  bloom.  Our  observations 
(Figs.  2  and  3)  do  not  support  this  hypothesis.  Further,  corre¬ 
lation  between  climatological  winds  and  ag  (0.63,  n  =  52) 
is  lower  than  between  wind  and  aph  (0.72,  n  =  52),  again 
suggesting  that  mixing  does  not  force  the  phase  shift. 

Are  the  lag  patterns  due  to  an  artifact  introduced  because  of 
the  constant  spectral  slope  S  =  0.02  nm-1  used  in  the  quasi- 
analytical  algorithm?  In  nature,  S  ranges  from  approximately 
0.012  nm-1  in  coastal  waters  to  0.034  nm-1  in  the  open 
ocean  (e.g.,  [1],  [9],  and  [25]).  Because  detrital  particles  tend 
to  have  similar  spectral  absorption  shapes  but  lower  spectral 
slope  values,  we  used  an  average  of  S  =  0.02  nm-1  in  our 
study,  which  is  . consistent  with  the  observations  made  in  the 
Sargasso  Sea  [25].  For  ag  <  0.01  m-1  (below  instrument 
detection  limits  for  443  nm)  S  can  only  be  inferred  from  the 
UV  wavelengths,  and  it  is  unclear  how  it  changes  when  CDOM 
is  photodegraded  (bleached).  Moran  et  al.  [24]  reported  that  for 
coastal  waters,  Suv  increased  by  6%  to  10%  when  CDOM  was 
60%  bleached.  We  therefore  performed  a  sensitivity  analysis 
by  changing  S  between  0.018  and  0.022  nm-1  either  randomly 
or  systematically  (i.e.,  increase  S  from  0.018  in  the  spring  to 
0.022  in  the  fall).  Results  show  that  the  magnitudes  of  aph  and 
ag  were  sensitive  to  the  choice  of  S,  but  the  temporal  patterns 
and  lags  between  these  properties  remained  independent  of 
these  changes  in  S. 

Carder  et  al.  [3]  had  documented  an  episode  of  ag  increase 
after  the  rapid  development  and  peaking  of  chlorophyll  and  aph 
in  a  coastal  phytoplankton  bloom  on  the  west  Florida  Shelf,  with 
a  lag  time  of  approximately  1.5  weeks.  Our  results  show  that 
similar  lags  exist  in  the  oligotrophic  ocean.  Loisel  et  al.  [17] 
showed  temporal  lag  patterns  in  the  particulate  backscattering 
coefficient  versus  chlorophyll  concentration  (band-ratio  algo¬ 
rithm)  from  SeaWiFS  in  the  North  Atlantic,  consistent  with  our 
findings  here.  Siegel  etal.  [31]  reported  global  CDM  patterns  as 
well  as  CDM  time-series  from  BATS.  Our  time-series  data  show 
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lower  values  but  similar  winter-spring  maxima  and  summer-fall 
minima  as  BATS.  Our  location  (~500  km  south  of  Bermuda) 
is  closer  to  the  gyre  center  and  is  characterized  by  clearer,  less 
productive  water.  Here,  the  uncertainties  in  the  magnitude  of 
aph  and  ag  estimates  are  likely  to  be  smaller,  which  helps  our 
study  of  the  relationship  between  phytoplankton  blooming  and 
CDOM  generation. 

Since  neither  vertical  mixing  nor  horizontal  advection,  nor  ar¬ 
tifacts  in  bio-optical  algorithms,  cause  the  lag  observed  between 
the  phytoplankton  and  CDOM,  we  present  the  following  hy¬ 
pothesis  to  explain  the  temporal  patterns  observed.  First,  CDOM 
is  produced  as  a  result  of  phytoplankton  grazing  and  cell  degra¬ 
dation,  causing  a  lag  in  the  a%  maxima.  Then,  CDOM  is  slowly 
destroyed  through  photodegradation  near  the  surface  [23],  As 
summer  starts,  phytoplankton  concentration  decreases  due  to 
nutrient  limitation  but  continues  to  generate  CDOM,  which  in 
turn  undergoes  bleaching  and  photodegradation  as  UV  light  in¬ 
tensity  increases  in  a  seasonal  intensification.  The  complex  in¬ 
teraction  between  these  processes  leads  to  an  apparently  slower 
rate  of  decrease  in  CDOM  during  the  summer  than  the  rate  of 
increase  after  the  winter-spring  phytoplankton  bloom. 

Our  results  (Figs.  2-4)  have  several  significant  implications. 
First,  they  provide  indirect  evidence  that  a  significant  portion 
of  CDOM  in  the  surface  layer  of  the  open  ocean  originates 
from  phytoplankton.  Phytoplankton  excretion,  bacterial  pro¬ 
cessing  of  decomposing  organic  matter,  and  sloppy  feeding 
by  zooplankton  are  all  possible  sources  of  CDOM.  During 
the  more  sparsely  sampled  in  situ  time  series  at  BATS,  the 
depth  maxima  of  the  dissolved  organic  carbon  (DOC)  and  the 
spring  bloom  of  diatomaceous  phytoplankton  coincided  [20], 
such  that  cell  lysis,  uncoupled  from  the  heterotrophic  sources 
of  zooplankton  grazing  and  bacterial  remineralization,  was 
inferred  as  the  major  source  of  DOC.  In  the  northern  Sargasso 
Sea,  macrozooplankton  dominate  during  spring,  and  limited 
sampling  suggests  that  herbivore  biomass  has  a  month  time  lag 
following  phytoplankton  [18].  Sampling  at  a  faster  rate  from 
daily  satellite  observations,  our  phase  shift  of  weeks,  rather 
than  months,  instead  suggests  that  marine  CDOM  may  mainly 
be  derived  from  lysis  of  phytoplankton  cells,  with  a  fraction 
derived  from  zooplankton  grazing  and  excretion. 

The  second  implication  is  the  impact  of  the  phase  shift  on  the 
bio-optical  band-ratio  algorithm  used  to  estimate  chlorophyll 
concentration  from  remote  sensing  observations  (for  oligotrophic 
waters  the  ratio  is  typically  between  R[S  at  443  and  555  nm,  as 
for  example  with  the  OC4v4  algorithm  of  [28]).  Clearly,  any 
bias  in  the  chlorophyll  estimates  is  of  fundamental  relevance  to 
estimating  primary  production  (see  [2]  and  references  therein). 
The  ratios  of  ag443/aph443  and  Chl(OC4)/aph443  clearly 
show  seasonal  variability  (Fig.  2).  This  implies  that  if  chlorophyll 
is  linearly  proportional  to  aph443,  then  the  band-ratio  OC4 
algorithm  will  overestimate  the  magnitude  and  possibly  the 
duration  of  the  spring  bloom  in  the  Sargasso  Sea.  We  estimate 
that  this  overestimate  of  chlorophyll  is  upward  of  10%  during 
the  spring.  Similarly,  there  is  an  underestimate  of  chlorophyll, 
during  the  fall  of  up  to  10%. 

The  bio-optical  artifact  introduced  by  this  seasonality  may 
require  reexamining  previous  inferences  about  biogeochemical 
and  physical  processes  in  the  subtropical  gyres.  For  example, 
color  anomalies  in  the  subtropical  gyres  have  been  hypothesized 
to  be  related  to  planetary  (Rossby)  waves.  Some  argue  that 


the  color  anomaly  may  be  due  to  local  phytoplankton  growth 
stimulated  by  wave-pumped  deep-waternutrients  [33],  but  others 
argue  that  the  cause  is  detrital  particles  advected  horizontally 
[5],  [15].  The  ability  to  separate  CDOM  from  chlorophyll,  the 
lag  observed  between  the  two,  and  the  new  knowledge  that  the 
lag  confounds  the  timing  and  magnitude  of  the  variability  in 
the  historical  OC4-derived  chlorophyll  concentration  suggests 
that  these  previous  hypotheses  require  further  research. 

Finally,  from  the  viewpoint  of  biogeochemical  modeling 
(e.g.,  [34]),  the  phase  shift  of  a%  relative  to  uph  for  the  deep 
sea  has  important  implications.  Deuser  et  al.  [6]  found  that 
seasonal  variability  in  surface  chlorophyll  concentration,  as 
derived  using  band  ratios  from  Coastal  Zone  Color  Scanner 
data,  was  reflected  in  the  sediment  flux,  and  that  flux  at  depth 
lagged  surface  production  by  about  one  and  a  half  months. 
If,  however,  prior  satellite  estimates  of  chlorophyll  were  not  a 
measure  of  the  actual  chlorophyll  concentration,  but  reflect  a 
convolution  of  lagged  pigment  and  CDOM  stocks,  the  conclu¬ 
sions  of  Deuser  et  al.  [6]  may  also  need  to  be  revised.  The  true 
time  difference  between  particle  export  and  surface  chlorophyll 
may  be  as  long  as  two,  or  two  and  one-half  months,  relative 
to  previous  estimates.  Particle  settling  rates  in  the  subtropical 
ocean  may  be  much  slower,  or  of  order  70  m  per  day,  instead  of 
the  100  m  per  day,  estimated  from  those  seminal  studies. 

IV.  Conclusion 

With  daily  high-resolution  SeaWiFS  data  and  a  recently 
developed  algorithm,  we  document  a  lag  (phase  shift)  between  the 
temporal  patterns  of  phytoplankton  pigments  and  CDOM  in  the 
North  Atlantic  subtropical  gyre.  Lag  times  are  about  two  weeks 
and  4-5  weeks  for  the  winter-spring  maxima  and  summer-fall 
minima,  respectively,  in  these  variables.  The  phase  shift  suggests 
that  CDOM  is  related  to  phytoplankton  growth,  grazing,  and 
degradation,  and  our  concept  of  linear  covariation  between  the 
two(thebasisfortheband-ratiochlorophyllalgorithms)mayneed 
to  be  revisited.  The  discovery  of  the  phase  shift  demonstrates  the 
significance  of  satellite  observations  in  process-related  studies. 
Advancing  and  using  this  technology  will  continue  to  move 
our  understanding  of  global  biogeochemical  cycles  forward. 
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